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Abstract. Frameshift mutations in protein-coding DNA sequences pro- 
duce a drastic change in the resulting protein sequence, which prevents 
classic protein alignment methods from revealing the proteins' common 
origin. Moreover, when a large number of substitutions are addition- 
ally involved in the divergence, the homology detection becomes difficult 
even at the DNA level. To cope with this situation, we propose a novel 
method to infer distant homology relations of two proteins, that accounts 
for frameshift and point mutations that may have affected the coding se- 
quences. We design a dynamic programming alignment algorithm over 
memory-efficient graph representations of the complete set of putative 
DNA sequences of each protein, with the goal of determining the two 
putative DNA sequences which have the best scoring alignment under a 
powerful scoring system designed to reflect the most probable evolution- 
ary process. This allows us to uncover evolutionary information that is 
not captured by traditional alignment methods, which is confirmed by 
biologically significant examples. 



1 Introduction 

In protein-coding DNA sequences, frameshift mutations (insertions or deletions 
of one or more bases) can alter the translation reading frame, affecting all the 
amino acids encoded from that point forward. Thus, frameshifts produce a dras- 
tic change in the resulting protein sequence, preventing any similarity to be 
visible at the amino acid level. 

When the coding DNA sequence is relatively well conserved, the similarity 
remains detectable at the DNA level, by DNA sequence alignment, as reported 
in several papers, including [1-4]. 

However, the divergence often involves additional base substitutions. It has 
been shown [5-7] that, in coding DNA, there is a base compositional bias among 
codon positions, that does not apply when the translation reading frame is 
changed. Hence, after a reading frame change, a coding sequence is likely to un- 
dergo base substitutions leading to a composition that complies with this bias. 
Amongst these substitutions, synonymous mutations (usually occurring on the 
third position of the codon) are more likely to be accepted by natural selection, 
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since they are silent with respect to the gene's product. If, in a long evolutionary 
time, a large number of codons in one or both sequences are affected by these 
changes, the sequence may be altered to such an extent that the common origin 
becomes difficult to observe by direct DNA comparison. 

In this paper, we address the problem of finding distant protein homologies, in 
particular when the primary cause of the divergence is a frameshift. We achieve 
this by computing the best alignment of DNA sequences that encode the target 
proteins. This approach relies on the idea that synonymous mutations cause 
mismatches in the DNA alignments that can be avoided when all the sequences 
with the same translation are explored, instead of just the known coding DNA 
sequences. This allows the algorithm to search for an alignment by dealing only 
with non-synonymous mutations and gaps. 

We designed and implemented an efficient method for aligning putative cod- 
ing DNA sequences, which builds expressive alignments between hypothetical 
nucleotide sequences that can provide some information about the common an- 
cestral sequence, if such a sequence exists. We perform the analysis on memory- 
efficient graph representations of the complete set of putative DNA sequences 
for each protein, described in Section 3.1. The proposed method, presented in 
Section 3.2, consists of a dynamic programming alignment algorithm that com- 
putes the two putative DNA sequences that have the best scoring alignment 
under an appropriate scoring system (Section 3.3) designed to reflect the actual 
evolution process from a codon-oriented perspective. 

While the idea of finding protein relations by frameshifted DNA alignments 
is not entirely new, as we will show in Section 2 in a brief related work overview, 
Section 4 - presenting tests performed on artificial data - demonstrates the 
efficiency of our scoring system for distant sequences. Furthermore, we validate 
our method on several pairs of sequences known to be encoded by overlapping 
genes, and on some published examples of frameshifts resulting in functional 
proteins. We briefly present these experiments in Section 5, along with a study of 
a protein family whose members present high dissimilarity on a certain interval. 
The paper is concluded in Section 6. 

2 Related Work 

The idea of using knowledge about coding DNA when aligning amino acid se- 
quences has been explored in several papers. 

A non- statistical approach for analyzing the homology and the "genetic semi- 
homology" in protein sequences was presented in [8,9]. Instead of using a sta- 
tistically computed scoring matrix, amino acid similarities are scored according 
to the complexity of the substitution process at the DNA level, depending on 
the number and type (transition/transversion) of nucleotide changes that are 
necessary for replacing one amino acid by the other. This ensures a differenti- 
ated treatment of amino acid substitutions at different positions of the protein 
sequence, thus avoiding possible rough approximations resulting from scoring 



them equally, based on a classic scoring matrix. The main drawback of this 
approach is that it was not designed to cope with frameshift mutations.. 

Regarding frameshift mutation discovery, many studies [1-4] preferred the 
plain BLAST [10,11] alignment approach: BLASTN on DNA and mRNA, or 
BLASTX on mRNA and proteins, applicable only when the DNA sequences are 
sufficiently similar. BLASTX programs, although capable of insightful results 
thanks to the six frame translations, have the limitation of not being able to 
transparently manage frameshifts that occur inside the sequence, for example by 
reconstructing an alignment from pieces obtained on different reading frames. 

An interesting approach for handling frameshifts at the protein level was 
developed in [12]. Several substitution matrices were designed for aligning amino 
acids encoded on different reading frames, based on nucleotide pair matches 
between respective codons. This idea has the advantage of being easy to use with 
any classic protein alignment tool. However, it lacks flexibility in gap positioning. 

On the subject of aligning coding DNA in presence of frameshift errors, some 
related ideas were presented in [13, 14] . The author proposed to search for protein 
homologies by aligning their sequence graphs (data structures similar to the ones 
we describe in Section 3.1). The algorithm tries to align pairs of codons, possibly 
incomplete since gaps of size 1 or 2 can be inserted at arbitrary positions. The 
score for aligning two such codons is computed as the maximum substitution 
score of two amino acids that can be obtained by translating them. This results 
in a complex, time costly dynamic programming method that basically explores 
all the possible translations. In Section 3.2, we present an algorithm addressing 
the same problem, more efficient since it aligns symbols, not codons, and more 
flexible with respect to scoring functions. Additionally, we propose to use a 
scoring system relying on codon evolution rather than amino acid translations, 
since we believe that, in frameshift mutation scenarios, the information provided 
by DNA sequence dynamics is more relevant than amino acid similarities. 

3 Our approach to distant protein relation discovery 

The problem of inferring homologies between distantly related proteins, whose 
divergence is the result of frameshifts and point mutations, is approached in this 
paper by determining the best pairwise alignment between two DNA sequences 
that encode the proteins. 

Given two proteins Pa and Pg, the objective is to find a pair of DNA se- 
quences, Da and Dg, such that translation{D a) — Pa and translation(D b) = 
Pb, which produce the best pairwise alignment under a given scoring system. 

The alignment algorithm (described in Section 3.2) incorporates a gap penalty 
that limits the number of frameshifts allowed in an alignment, to comply with the 
observed frequency of frameshifts in a coding sequence's evolution. The scoring 
system (Section 3.3) is based on possible mutational patterns of the sequences. 
This leads to reducing the false positive rate and focusing on alignments that 
are more likely to be biologically significant. 



3.1 Data structures 



An explicit enumeration and pairwise alignment of all the putative DNA se- 
quences is not an option, since their number increases exponentially with the 
protein's length 1 . Therefore, we represent the protein's "back-translation" (set 
of possible source DNAs) as a directed acyclic graph, whose size depends linearly 
on the length of the protein, and where a path represents one putative sequence. 

As illustrated in Figure 1(a), the graph is organized as a sequence of length 
3n where n is the length of the protein sequence. At each position i in the 
graph, there is a group of nodes, each representing a possible nucleotide that 
can appear at position i in at least one of the putative coding sequences. Two 
nodes at consecutive positions are linked by arcs if and only if they are either 
consecutive nucleotides of the same codon, or they are respectively the third and 
the first base of two consecutive codons. No other arcs exist in the graph. 

Note that in the implementation, the number of nodes is reduced by using 
the IUPAC nucleotide codes. If the amino acids composing a protein sequence 
are non-ambiguous, only 4 extra nucleotide symbols - R,Y , H and N — are nec- 
essary for their back-translation. In this condensed representation, the number 
of ramifications in the graph is substantially reduced, as illustrated by Figure 1. 
More precisely, the only amino acids with ramifications in their back-translation 
are amino acids R, L and S, each encoded by 6 codons with different prefixes. 

3.2 Alignment algorithm 

We use a dynamic programming method, similar to the Smith- Waterman algo- 
rithm, extended to data structures described in Section 3.1 and equipped with 
gap related restrictions. 

Given the input graphs Ga and Gb obtained by back-translating proteins Pa 
and Pb, the algorithm finds the best scoring local alignment between two DNA 
sequences comprised in the back-translation graphs (illustrated in Figure 2). 
The alignment is built by filling each entry M[i,j, («A,as)] of a dynamic pro- 
gramming matrix M, where i and j are positions of the first and second graph 
respectively, and (a a, as) is a pair of nodes that can be found in Ga at position 
i, and in Gb at position j, respectively. An example is given in Figure 3. 

The dynamic programming algorithm begins with a classic local alignment 
initialization (0 at the top and left borders), followed by the recursion step de- 
scribed in equation (1). The partial alignment score from each cell M[i,j, (ct/i, «b)] 
is computed as the maximum of 6 types of values: 

(a) (similarly to the classic Smith- Waterman algorithm, only non-negative 
scores are considered for local alignments). 

(b) the substitution score of symbols («A,as), denoted score{aA 1 as) 1 added 
to the score of the best partial alignment ending in M[i — l,j — 1], provided 
that the partially aligned paths contain a a on position i and as on position 

1 With the exception of M and W, which have a single corresponding codon, all amino 
acids are encoded by 2, 3, 4 or 6 codons. 



Fig. 1. Example of fully repre- 
sented (a) and condensed (b) 
back-translation graph for the 
amino acid sequence YSH. 



Fig. 2. Alignment example. A path (correspond- 
ing to a putative DNA sequence) was chosen from 
each graph so that the match/mismatch ratio is 
maximized. 
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■ M[i, j] is a "ceZZ" of M corresponding to 
position i of the first graph and position 
j of the second graph. 
M[i,j] contains entries (aA,as) corre- 
sponding to pairs of nodes occurring in 
the first graph at position i, and in the 
second graph at position j , respectively. 



Fig. 3. Example of dynamic programming matrix M. 



j respectively; this condition is ensured by restricting the entries of M[i — 
1, j — 1] to those labeled with symbols that precede a a and as in the graphs. 

(c) the cost singleGap Penalty of a frameshift (gap of size 1 or extension of a 
gap of size 1) in the first sequence, added to the score of the best partial 
alignment that ends in a cell M[i, j — 1, (a a, /3b)], provided that /3g precedes 
as in the second graph; this case is considered only if the number of allowed 
framcshifts on the current path is not exceeded, or a gap of size 1 is extended. 

(d) the cost of a frameshift in the second sequence, added to a partial alignment 
score defined as above. 

(e) the cost tripleGap Penalty of removing an entire codon from the first se- 
quence, added to the score of the best partial alignment ending in a cell 

M[i,j-3,(a A ,pB)]. 

(f ) the cost of removing an entire codon from the second sequence, added to the 
score of the best partial alignment ending in a cell M[i — 3, j, {Pa, c*b)] 

We adopted a non-monotonic gap penalty function, which favors insertions 
and deletions of full codons, and does not allow a large number of frameshifts 
- very rare events, usually eliminated by natural selection. As can be seen in 
equation (1), two particular kinds of gaps are considered: i) frameshifts - gaps 
of size 1 or 2, with high penalty, whose number in a local alignment can be 



limited, and ii) codon skips - gaps of size 3 which correspond to the insertion 
or deletion of a whole codon. 



M[i,j, (a A ,a B )] = 

'0 (a) 

M[i - 1, j - 1, (/3a, Pb)] + score(a A ,a B ), Pk € pred(a k ); (b) 

max < ~ X ' ( a ^' + singleGapPenalty) , /3 B G pred(a B ); (c) (1) 

j (M[i — 1, j, (/Sa^b)] + singleGapPenalty) , (3 A G prec^a^); (d) 

(M[i,j — 3, (o^jAb)] + tripleGapPenalty) , j > 3 (e) 

k (M[i - 3, j, (13 a, a B )] + tripleGapPenalty) , i > 3 (f) 



3.3 Translation-dependent scoring function 

In this section, we present a new translation-dependent scoring system suitable 
for our alignment algorithm. The scoring scheme we designed incorporates in- 
formation about possible mutational patterns for coding sequences, based on 
a codon substitution model, with the aim of filtering out alignments between 
sequences that are unlikely to have common origins. 

Mutation rates have been shown to vary within genomes, under the influence 
of several factors, including neighbor bases [15]. Consequently, a model where 
all base mismatches are equally penalized is oversimplified, and ignores possibly 
precious information about the context of the substitution. 

With the aim of retracing the sequence's evolution and revealing which base 
substitutions are more likely to occur within a given codon, our scoring system 
targets pairs of triplets (a,p, a), were a is a nucleotide, p is its position in the 
codon, and a is the amino acid encoded by that codon, thus differentiating 
various contexts of a substitution. There are 99 valid triplets out of the total of 
240 hypothetical combinations. 

Pairwise alignment scores are computed for all possible pairs of valid triplets 
(*i,*2) = ((oiiiPij ai), (ct2,P2, 0,2)) as a classic log-odds ratio: 

score(t 1 ,t 2 ) = A log (2) 

0*1*2 

where f tl t 2 is the frequency of the t\ <-> t 2 substitution in related sequences, and 
k*i*2 = p(ti)p(t2) is the background probability. 

In order to obtain the foreground probabilities /t^, we will consider the 
following scenario: two proteins are encoded on the same DNA sequence, on 
different reading frames; at some point, the sequence was duplicated and the 
two copies diverged independently; we assume that the two coding sequences 
undergo, in their independent evolution, synonymous and non-synonymous point 
mutations, or full codon insertions and removals. 

The insignificant amount of available real data that fits our hypothesis does 
not allow classical, statistical computation of the foreground and background 
probabilities. Therefore, instead of doing statistics on real data directly, we will 
rely on codon frequency tables and codon substitution models. 



We assume that codon substitutions in our scenarios can be modeled by 
a Markov model presented in [16] 2 which specifies the relative instantaneous 
substitution rate from codon i to codon j as: 

if i or j is a stop codon, or 

if z — >- j requires more than 1 nucleotide substitution, 
n I TTj if i — > j is a synonymous transversion, 
3 | TTjK if i — > j is a synonymous transition, 

irjU> if i — > j is a nonsynonymous transversion, 
^ nj kuj if i — > j is a nonsynonymous transition. 

for all i ^ j. Here, the parameter oj represents the nonsynonymous-synonymous 
rate ratio, k the transition-transversion rate ratio, and ttj the equilibrium fre- 
quency of codon j. As in all Markov models of sequence evolution, absolute rates 
are found by normalizing the relative rates to a mean rate of 1 at equilibrium, 
that is, by enforcing ^ i Sj^i ^iQij = 1 an d completing the instantaneous rate 
matrix Q by defining Qu — — J2j^i Qij to & ve a form in which the transition 
probability matrix is calculated as P{9) — e e ® [18]. Evolutionary times 9 are 
measured in expected number of nucleotide substitutions per codon. 

With this codon substitution model, j ti tj can be deduced in several steps. 
Basically, we first need to identify all pairs of codons with a common subse- 
quence, that have a perfect semi-global alignment (for instance, codons CAT 
and ATG satisfy this condition, having the common subsequence AT; this ex- 
ample is further explained below). We then assume that the codons from each 
pair undergo independent evolution, according to the codon substitution model. 
For the resulting codons, we compute, based on all possible original codon pairs, 
p((ai,pi, q), (aj,pj, cj)) - the probability that nucleotide a^, situated on posi- 
tion pi of codon a , and nucleotide ctj , situated on position pj of codon Cj have 
a common origin (equation (5)). From these, we can immediately compute, as 
shown by equation (6), p((ai,Pi,di), (aj,pj,aj)), corresponding in fact to the 

foreground probabilities /t^, where ti = (m.pi.ai) and tj — (aj 7 pj 7 aj). 



In the following, p(ci — > C2) stands for the probability of the event codon c\ 
mutates into codon C2 in the evolutionary time 0, and is given by P Cl , C2 (8). 

Ci [interval].] = c 2 [interval 2 ] states that codon ci restricted to the positions 
given by intervali is a sequence identical to ci restricted to interval- This 
is equivalent to having a word w obtained by "merging" the two codons. For 
instance, if c\ = CAT and c 2 = ATG, with their common substring being 
placed in intervali = [2. .3] and interval = [1-.2] respectively, w is CATG. 

Finally, p(ci [intervali] = c 2 [interval 2 ]) is the probability to have c\ and c 2 , 
in the relation described above, which we compute as the probability of the word 
w obtained by "merging" the two codons. This function should be symmetric, it 

2 Another, more advanced codon substitution model, targeting sequences with over- 
lapping reading frames, is proposed and discussed in [17]. It does not fit our scenario, 
because it is designed for overlapping reading frames, where a mutation affects both 
translated sequences, while in our case the sequences become at one point indepen- 
dent and undergo mutations independently. 



(3) 



should depend on the codon distribution, and the probabilities of all the words 
to of a given length should sum to 1. However, since we consider the case where 
the same DNA sequence is translated on two different reading frames, one of the 
two translated sequences would have an atypical composition. Consequently, the 
probability of a word w is computed as if the sequence had the known codon 
composition when translated on the reading frame imposed by the first codon, 
or on the one imposed by the second. This hypothesis can be formalized as: 

p(w) — p(w on rfi OR w on r/2) = p r ^ (w) + p 1 ^ 2 (w) — p 7 "* 1 (w) ■ p r ^ 2 (w) (4) 

where p 7 "' 1 (w) and p r -> 2 (w) are the probabilities of the word w in the reading 
frame imposed by the position of the first and second codon, respectively. This is 
computed as the products of the probabilities of the codons and codon pieces that 
compose the word w in the established reading frame. In the previous example, 
the probabilities of w = CATG in the first and second reading frame are: 

p rfl (CATG) = p'CAT) ■ p(G * *) = p(CAT) ■ ^ p(c) 

c:c starts with G 

p rh (CATG) = p(* * C) ■ p(ATG) = p( c ) ' p( at G) 

c:c ends with C 

The values of p((ai 7 pi, a), (ctj 7 pj 7 Cj)) are computed as: 

c'j[intervalj]) ■ p(c[ A q) ■ p(c'j A cj) (5) 



from which obtaining the foreground probabilities is straightforward: 

futj =p((ai,Pi,a i ),(a j ,p j ,a j )) = p((ai,p u c,), (aj,pj, c,)) (6) 

c^encodes ai , 
cj encodes aj 

The background probabilities of (ti,tj), b^u, can be simply expressed as 
the probability of the two symbols appearing independently in the sequences: 

but, = &(a i ,Pi,o 4 ),(ai,Pi,a J ) = ^ 

□^encodes a i , 
c^encodes aj 

Substitution matrix for ambiguous symbols From matrices built as explained 
above, the versions that use IUPAC ambiguity codes for nucleotides (as pro- 
posed in the final paragraph of 3.1) can be computed: the score of pairing two 
ambiguous symbols is the maximum over all substitution scores for all pairs of 
nucleotides from the respective sets. 

Score evaluation The score significance is estimated according to the Gumbel 
distribution, where the parameters A and K arc computed with the method 
described in [19,20]. Since the forward alignment and the reverse complemen- 
tary alignment are two independent cases with different score distributions, two 
parameter pairs, Xf w ,Kf w and A rc , K rc are computed and used in practice. 



p(c\ [intervali 



.'..c'j :a' z [i„ter»al 4 ] = cj [inter,.!, ] 

Pi£intervali : pj ^interval j 



4 Validation 



To validate the translation-dependent scoring system we designed in the previous 
section, we tested it on an artificial data set consisting in 96 pairs of protein se- 
quences of average length 300. Each pair was obtained by translating a randomly 
generated DNA sequence on two different reading frames. Both sequences in each 
pair were then mutated independently, according to codon mutation probability 
matrices corresponding to each of the evolutionary times 0.01, 0.1, 0.3, 0.5, 0.7, 
1.0, 1.5, 2.00 (measured in average number of mutations per codon). 

To this data set we applied four variants of alignment algorithms: i) classic 
alignment of DNA sequences using classic base substitution scores and afiinc gap 
penalties; ii) classic alignment of DNA sequences using a translation-dependent 
scoring scheme designed in Section 3.3; hi) alignment of back-translation graphs 
(Section 3.2) using classic base substitution scores and afhne gap penalties; 
iv) alignment of back-translation graphs using a translation-dependent scoring 
scheme. For the tests involving translation-dependent scores, we used scoring 
functions corresponding to evolutionary times from 0.30 to 1.00. 

Table 1 briefly shows the e- values of the scores obtained with each setup when 
aligning sequence pairs with various evolutionary distances. While all variants 
perform well for highly similar sequences, we can clearly deduce the ability of the 
translation-dependent scores to help the algorithm build significant alignments 
between sequences that underwent important changes. 
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Table 1. Order of the e- values of the scores obtained by aligning artificially diverged 
pairs of proteins resulted from the translation of the same ancestral sequence on two 
reading frames. W TDS < evolutionary distance> = translation-dependent scores; clas- 
sic substitution scores: match = 3, transversion = -4, transition = -2. 



The resulting alignments reveal that, even after many mutations, the transla- 
tion-dependent scores manage to recover large parts of the original shared se- 
quence, by correctly aligning most positions. On the other hand, with classic 
match/mismatch scores, the algorithm usually fails to find these common zones. 
Moreover, due to the large number of mismatches, the alignment has a low score, 



comparable to scores that can be obtained for randomly chosen sequences. This 
makes it difficult to establish whether the alignment is biologically meaningful 
or it was obtained by chance. This issue is solved by the translation-dependent 
scores by uneven substitution penalties, according to the codon mutation models. 

We conclude that the usage of translation-dependent scores makes the algo- 
rithm more robust, able to detect the common origins even after the sequences 
underwent many modifications, and also able to filter out alignments where the 
nucleotide pairs match by pure chance and not due to evolutionary relations. 

5 Experimental results 

5.1 Tests on known overlapping and frameshifted genes 

We tested the method on pairs of proteins known to be encoded by overlap- 
ping genes in viral genomes (phage X174 and Influenza A) and in E.coli plas- 
mids, as well as on the newly identified overlapping genes yaaW and htgA from 
E.coli K12 [21]. In all cases, we obtained perfect identification of gene overlaps 
with simple substitution scores and with translation-dependent scoring matrices 
corresponding to low evolutionary distances (at most 1 mutation per codon ). 
Translation-dependent scoring matrices of higher evolutionary distances favor, 
in some (rare) cases, substitutions instead of matches within the alignment. This 
is a natural consequence of increasing the codon's chance to mutate, and it il- 
lustrates the importance of choosing a score matrix corresponding to the real 
evolutionary distance. Our method was also able to detect, directly on the pro- 
tein sequences, the frameshifts resulting in functional proteins reported in [1-4]. 

5.2 New divergence scenarios for orthologous proteins 

In this section we discuss the application of our method to FMR1NB (Fragile 
X mental retardation 1 neighbor protein) family. The Ensembl database [22] 
provides 23 members of this family, from mammalian species, including human, 
mouse, dog and cow. Their multiple alignment, provided by Ensembl, shows 
high dissimilarity on the first part (100 amino acids approximately), and good 
conservation on the rest of the sequence. We apply our alignment algorithm on 
proteins from several organisms, where the complete sequence is available. 

We performed our experiments with translation-dependent scoring matrices 
corresponding to 0.3, 0.5 and 0.7 mutations per codon. Given that, in our sce- 
nario (presented in section 3.3), the divergence is applied on two reading frames, 
this implies an overall mutation rate of 0.6, 1.0 and 1.4 mutations per codon re- 
spectively. Thus, the mutation rate per base reflected by our scores is less than 
0.5, which is approximately the nucleotide substitution rate for mouse relative 
to human [23]. The number of allowed frameshifts was limited to 3. The gap 
penalties were set in all cases to -20 for codon indels, -20 for size 1 gaps and -5 
for the extension of size 1 gaps (size 1 and size 2 gaps correspond to frameshifts). 
These choices were made so that the penalty for codon indels is higher than the 
average penalty for 3 substitutions. 
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Fig. 4. Human and mouse FMR1NB proteins, aligned using a translation-dependent 
matrix of evolutionary distance 0.7 (the sign of each substitution score appears on the 
fourth row). The size 4 gap corresponds to a frameshift that corrects the reading frame. 

Figure 4 presents a fragment of the alignment obtained on the FMR1NB 
proteins of human (gene ID ENSG00000176988) and mouse (gene ID ENS- 
MUSG00000062170). The algorithm finds a frameshift near the 100th amino 
acid, managing to align the initial part of the proteins at the DNA level. Similar 
frameshifted alignments are obtained for human vs. cow and human vs. dog, 
while alignments between proteins of primates do not contain frameshifts. The 
consistency of the frameshift position in these alignments supports the evidence 
of a frameshift event that might have occurred in the primate lineage. 

If confirmed, this frameshift would have modified the first topological domain 
and the first transmembrane domain of the product protein. Interestingly, the 
FMR1NB gene occurs nearby the Fragile X mental retardation 1 gene (FMR1), 
involved in the corresponding genetic disease [24] . 

6 Conclusions 

In this paper, we addressed the problem of finding distant protein homologies, in 
particular affected by frameshift events, from a codon evolution perspective. We 
search for protein common origins by implicitly aligning all their putative coding 
DNA sequences, stored in efficient data structures called back-translation graphs. 
Our approach relies on a dynamic programming alignment algorithm for these 
graphs, which involves a non-monotonic gap penalty that handles differently 
frameshifts and full codon indels. We designed a powerful translation-dependent 
scoring function for nucleotide pairs, based on codon substitution models, whose 
purpose is to reflect the expected dynamics of coding DNA sequences. 

The method was shown to perform better than classic alignment on artificial 
data, obtained by mutating independently, according to a codon substitution 
model coding sequences translated with a frameshift. Moreover, it successfully 
detected published frameshift mutation cases resulting in functional proteins. 

We then described an experiment involving homologous mammalian proteins 
that showed little conservation at the amino acid level on a large region, and 
provided possible frameshifted alignments obtained with our method, that may 
explain the divergence. As illustrated by this example, the proposed method 
should allow to better explain a high divergence of homologous proteins and to 
help to establish new homology relations between genes with unknown origins. 

An implementation of our method is available at http://bioinfo.lifl.fr/path/. 
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